Monogamy inequality for distributed Gaussian entanglement 
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We show that for all zi-mode Gaussian states of continuous variable systems, the entanglement shared among 
n parties exhibits the fundamental monogamy property. The monogamy inequality is proven by introducing 
the Gaussian tangle, an entanglement monotone under Gaussian local operations and classical communication, 
which is defined in terms of the squared negativity in complete analogy with the case of rc-qubit systems. Our 
results elucidate the structure of quantum correlations in many-body harmonic lattice systems. 
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Quantum entanglement, at the heart of quantum correla- 
tions and a direct consequence of the superposition principle, 
cannot be freely shared among many parties, unlike classi- 
cal correlations. This is the so-called monogamy property [1] 
and is one of the fundamental traits of entanglement and of 
quantum mechanics itself [2]. Seminal observations on the 
monogamy property and its precise quantitative statement in 
mathematical terms are due to Coffman, Kundu, and Wootters 
^ . They proved the following inequality for a generic state 
Pabc of three qubits, 



T(PA:BC) > t(PA; B ) + T(p A :c) , 



(1) 



where t is an entanglement monotone |01 known as the tangle 
JH, t(pa.bc) stands for the bipartite entanglement across the 
bipartition A : BC, and Pa-.b(q = ^c(B){Pabc)- Inequality (HJ 
clearly elucidates the restriction on the sharing of entangle- 
ment among the three parties. The monogamy inequality ([TJ 
holds for the tangle defined through the square of the concur- 
rence Hi. Recently, Osborne and Verstraete [7] have gen- 
eralized the monogamy inequality to n-qubit systems, prov- 
ing a longstanding conjecture formulated in Ref. y|], with im- 
portant consequences for the description of the entanglement 
structure in many-body spin systems. 

For higher-dimensional systems much less is known on the 
qualification, let alone the quantification of bipartite and mul- 
tipartite entanglement, the situation worsening with increas- 
ing dimension of the Hilbert space due to the exponential in- 
crease in the complexity of states. Remarkably, in the limit of 
continuous variable (CV) systems with an infinite dimensional 
Hilbert space, if one focuses on the theoretically and practi- 
cally relevant class of Gaussian states an almost as compre- 
hensive characterization of entanglement has been achieved 
as in the case of qubit systems [8]. In this context, the natu- 
ral question arises whether the monogamy inequality holds as 
well for entanglement sharing in CV systems, and in particular 
in generic Gaussian states of harmonic lattices. The first step 
towards answering this question has been taken by Adesso and 
Illuminati [9] (see also lldO . They proved the monogamy 



inequality for arbitrary three-mode Gaussian states and for 
symmetric n-mode Gaussian states, defining the CV tangle or 
"contangle" as the square of the logarithmic negativity (an en- 
tanglement monotone |l ill ). However, it is known [12] that in 
two-qubit systems the concurrence is equivalent to another re- 
lated entanglement measure, the negativity lTl3tl . so it appears 
natural to promote the tangle to CV systems by defining it in 
terms of the squared negativity itself 11411 . 

In this paper we provide the complete answer to the ques- 
tion posed above. We prove that the monogamy inequality 
does hold for all Gaussian states of multimode CV systems 
with an arbitrary number n of modes and parties A\,. . . ,A„, 
thus generalizing the results of |9l llOfl. As a measure of bi- 
partite entanglement, we define the Gaussian tangle via the 
square of negativity, in direct analogy with the case of n-qubit 
systems [7]. Our proof is based on the symplectic analysis of 
covariance matrices and on the properties of Gaussian mea- 
sures of entanglement The monogamy constraint has 
important implications on the structural characterization of 
entanglement sharing in CV systems |0, [hJ, in the context 
of entanglement frustration in harmonic lattices flitl . and for 
practical applications such as secure key distribution and com- 
munication networks with continuous variables. 

For a A] : A% . . .A„ bipartition associated to a pure Gaus- 
sian state p ( ^ B with A = A\ (a subsystem of a single mode) 
and B — A2 . . . A„, we define the following quantity 



(2) 



Here, N(p) = (||p r ' 1 || 1 - l)/2 is the negativity fill [1311, 
denotes the trace norm, and p TA stands for the partial transpo- 
sition of p with respect to the subsystem A. The functional tq, 
like the negativity N, vanishes on separable states and does 
not increase under local operations and classical communica- 
tion (LOCC), i.e., it is a proper measure of pure-state bipartite 
entanglement |4[]. It can be naturally extended to mixed Gaus- 
sian states pa.b via the convex roof construction 



TciPA-.B) = inf }.PiTGipf\ 



(3) 



2 



where the infimum is taken over all convex decompositions of 
Pa.b in terms of pure Gaussian states p^ p) : pa.b = H; PiP- ■ By 
virtue of the convex roof construction, tq [Eq. ©J is an entan- 
glement monotone under Gaussian LOCC (GLOCC) Jl5 , 171. 

Henceforth, given an arbitrary n-mode Gaussian state 
PA!-a 2 ..a„> we refer to tq [Eq. 0)] as the Gaussian tangle and 
we now prove the general monogamy inequality 



where a = ( r " y ' 2 ) is the CM of the single-mode reduced 
Gaussian state p Al = T*A 2 ..A n pA i A 2 ..A n - The CM y of a n- 
mode pi<r<? Gaussian state is characterized by the condition 
II 1 711 -J n yJ n y = l2« > which implies 



det a + V det 5/ = 1 



1=2 



(7) 



/=2 



(4) 



To this end, we can assume without loss of generality that 
the reduced two-mode states p Al A, = Tu 2 ,.A,^Ai +i -A n PAiA 2 ..A„ 
of subsystems (A\Ai) (I = 2, ...,n) are all entangled. In 
fact, if for instance Pa x a 2 is separable, then tg(PAiA 3 ,.A„) ^ 
tg(PAiA 2 ..a„) because the partial trace over the subsystem A 2 
is a local Gaussian operation that does not increase the Gaus- 
sian entanglement. Furthermore, by the convex roof con- 
struction of the Gaussian tangle, it is sufficient to prove the 
monogamy inequality for any pure Gaussian state pj^. A A 
(see also Refs. 0, 01). Therefore, in the following we 
can always assume that Pa,Ai..a„ is a pure Gaussian state for 
which the reduced states pAyA, (I - 2, ... ,n) are all entangled. 

Some technical preliminaries are in order. A n-mode (pure 
or mixed) Gaussian state p is completely characterized by the 
covariance matrix (CM) y^ = 2Ti[p(Rj-dj)(Rk-dk)]-i(J»)jk, 
and by the displacement vector d% = TxipRk). Here R = 



, 1/2.-, -1/2 D 



1/2 —1/2 T 

Qn, w„ P n ) , with Qk and Pk the 



1 iii^j 

canonical quadrature-phase operators (position and momen- 
tum) for mode k with energy u>k, and J„ = ®"j =l J\ with 

J 1 - ( °j {jV the symplectic matrix fl. Every unitary transfor- 
mation for general (pure or mixed) Gaussian states p i-> U pU ' 
is associated to a symplectic transformation y h-> SyS T with 
S € Sp(2n,R) = {5|57„5 r = J n ). Positivity of the density 
matrix p is expressed in terms of y as 



y + z7„ > . 



(5) 



Under partial transposition, PayA 2 ..a„ p^-A a ' tne 

CM 7 is transformed to 7 = F„yF n , where F„ - 
diag(l,-l,l,l,...,l). 

We start by computing the left-hand side of Eq. ©. Since 
PAyAi..A n is a 1 X (n - 1) pure Gaussian state, it can be trans- 
formed as iflilfll 



AiA 2 ..A n U = pa, ® Pa;, ® • ■ 



5 a*;, 



(6) 



by a local unitary transformation U = Ua, ® Ua 2 ..a„ without 
changing the amount of entanglement across the bipartition 
A\ : A2 . . .A n . In the right-hand side of Eq. ©, Pa, : a; is a 
pure two-mode Gaussian state (a two-mode squeezed state) 
and pa; (I = 3, ... ,«) are vacuum states. Thus, t g (paiA^..a„) 
is equal to TGip Al A' 2 ) = N 2 (p Al A 2 )- In turn, N 2 (p Al A' 2 ) = 



(vj - 1) 2 /4 111 lLl20ll . where v_ denotes the smallest symplectic 
eigenvalue of y Al A' 2 = F 2 y Al A' 2 F 2 , with being the CM 

of Pai:a;- It is easy to compute v_; v_ = Vdeta - Vdeta - 1, 



where 61 is the matrix encoding intermodal correlations be- 
tween mode 1 and mode I in the reduced state Pa^a, (I = 
2, . . . , n), described by a CM 



7 A, A, 
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ri,2 


7l,2/-l 


7\,2l s 


T2,l 


T2,2 


72,2/- 1 


72,21 


72/- 1,1 


72/- 1,2 


72/-l,2/-l 


72/- 1,2/ 


* 721, 1 


72Z.2 


72/.2/-1 


72/,2/ ) 



Pi 



(8) 

As pAyA, is entangled, det 5/ is negative [21]. It is useful to 
introduce the auxiliary quantities A/ = -4 det 5/ > 0. The 
Gaussian tangle for Pa i a 2 ..a„ is then written as 



T G (PA I :A 2 ...A„)=^(V: 1 -1) 2 =/ 



where /(f) = (g-\f)- 1/2) 2 , with #(?) 



/ n 
V/=2 



(9) 



VFT4- Vf. (io) 



We observe that /(f) /f is an increasing function for f > and 
/(0) = so / is a star-shaped function: /(cf) < c/(f) for 
c e [0, 1] and f > Q. Therefore, we have /(f) < ^/(f + s) 
and/(i)< -f-f(t+s) for f, s > to obtain /(f) + f(s) < f(t+s). 
That is, / is superadditive l23ll . Hence |24], 



/ 



v 1=2 



2> *E*4>- 



(id 



/=2 



Each term in the right-hand side is well defined since A/ > 0. 

We are now left to compute the right-hand side of Eq. (01, 
i.e. the bipartite entanglement in the reduced (mixed) two- 
mode states Pa,:a, (I — 2, ...,n). We will show that the 
corresponding Gaussian tangle is bounded from above by 
/(A;), which will therefore prove the monogamy inequality 
via Eq. (fTTT) . To this aim, we recall that any bipartite and mul- 
tipartite entanglement in a Gaussian state is fully specified in 
terms of its CM, as the displacement vector of first moments 
can be always set to zero by local unitary operations, which 
preserve entanglement by definition. It is thus convenient to 
express the Gaussian tang le directly in terms of the CMs. Fol- 
lowing Refs. EES [Sill, the definition [Eq. ©] for the 
Gaussian tangle of a mixed Gaussian state with CM 7a, : a, can 
be rewritten as 



TO 



(yA 1 Aj ) = inf {r G (7^)lr^, ^ 7m a] , (12) 

Ta, a. 



where the infimum is taken over all CMs 7 A p ^. Al of pure Gaus 



sian states such that 7a, -a, > 7 Ai a • The quantities A/ and 



3 



Tc(yA,:A,) for any /, as well as every single-mode reduced 
determinant, are Sp(2, R) ffi " -invariants. For each two-mode 
partition described by Eq. (|8), we can exploit such local- 
unitary freedom to put the CM y A] - Al in standard form [25] 
with a = diag{a, a), Bi = diag{/>, b}, and 61 = diag{c + , c_], 
where c+ > \cJ\ 112 ll 12611 . The condition [Eq. (0] for yA^A, is 



thus equivalent to the following inequalities 

a>\,b>\,ab-c 2 ± >\; (13) 
Aely A[ A, + l=(ab- c 2 + ){ab - c 2 ) + 1 > a 2 + b 2 + 2c + c_.(14) 



Furthermore, since the state Pa,Ai is entangled, we have 1121 1 



(ab - c 2 + ){ab - c 2 ) + 1 < a 2 + b 2 - 2c + c_ . (15) 

From Eqs. ( fT4] > and JIBl . it follows that c_ < 0. In Eq. ( fT2l , 
t g(7 a Ia ) - f(4deta (p> - 4), which is an increasing function 
of deta (p \ where is the first 2x2 principal submatrix 
of yf^A,- The infimum of the right-hand side of Eq. (fT2l is 



achieved by the pure-state CM y A ^. A (with y A '. Al < yA^A, and 
"fAi-A + — ^) tnat minimizes deta^. The minimum value 



of det is given by mino<e<2 ff 01511 . where: m(ff) = 1 + 
h 2 (ff)/h 2 (ff) , with «i(6>) = £_ + yftjcosO, and « 2 (6») = 2(a/> - 

c 2 )(a 2 + b 2 + 2c + c_) - (ft/ ^) cos (9+ (a 2 - £ 2 ) a/1 ~£+/Y sin 0. 
Here 



ftt = c + (ab — c_) + C- , 
T} = [a — b(ab — c 2 )] [Z? — a(a/b - c 2 )] . 
£ = 2afrc 3 . + (a 2 + fr 2 )c+c 2 



(16) 
(17) 



+ [a 2 + b 2 - 2a 2 b 2 ]c- - ab(a 2 + b 2 - 2)c + . (18) 

Moreover, m(n) > mino<e<2 ff m{8) and therefore 

T G (y Al A,) < f(4m(n) - 4) = / (4 tf/fr) , (19) 

where £\ = h\{n) and t^i = h2{n). Finally, one can prove that 
(see the Appendix) 



A, = -4det<J/ = -4c + c_ > 4 f 2 /ft , 



(20) 



which, being /(f) [Eq. ( [TOl iI an increasing function of t, entails 
that /(A/) > /(4 £ 2 /£>). Combining this with Eq. (15) leads 
to the crucial Sp(2, R) e "-invariant condition 

T G (rA,:A,)</(A,), (21) 

which holds in general for all / = 2 . . . n and does not rely 
on the specific standard form of the reduced CMs y Al A, l25ll . 
Then, recalling Eqs. (O, ( fTTT i. and ( 1211 ). Inequality (|4]i is estab- 
lished. This completes the proof of the monogamy constraint 
on CV entanglement sharing for pure n-mode Gaussian states 
distributed among n parties. As already mentioned, the proof 
immediately extends to arbitrary mixed Gaussian states by the 
convexity of the Gaussian tangle [Eq. ©]. ■ 
Summarizing, we have defined the Gaussian tangle tq, an 
entanglement monotone under GLOCC, and proved that it is 
monogamous for all multimode Gaussian states distributed 



among multiple parties. The implications of our result are 
manifold. The monogamy constraints on entanglement shar- 
ing are essential for the security of C V quantum cryptographic 
schemes 12711 . because they limit the information that might 
be extracted from the secret key by a malicious eavesdropper. 
Monogamy is useful as well in investigating the range of cor- 
relations in Gaussian matrix-product states of harmonic rings 
I2H1 . and in understanding the entanglement frustration occur- 
ring in ground states of many-body harmonic lattice systems 
B16I1 . which, following our findings, may be now extended to 
arbitrary states beyond symmetry constraints. 

On the other hand, investigating the consequences of the 
monogamy property on the structure of entanglement sharing 
in generic Gaussian states along the lines of Refs. JSJ llOll . 
reveals that there exist states that maximize both the pair- 
wise entanglement in any reduced two-mode partition, and 
the residual distributed (multipartite) entanglement obtained 
as a difference between the left-hand and the right-hand side in 
Eq. (0). The simultaneous monogamy and promiscuity of CV 
entanglement (unparalleled in qubit systems), which can be 
unlimited in four-mode Gaussian states J29I1 . allows for novel, 
robust protocols for the processing and transmission of quan- 
tum and classical information II 1 CXI . The monogamy inequality 
[Eq. ©J bounds the persistency of entanglement when one or 
more nodes in a CV communication network sharing generic 
n-mode Gaussian resource states are traced out. 

At a fundamental level, the proof of the monogamy prop- 
erty for all Gaussian states paves the way to a proper quantifi- 
cation of genuine multipartite entanglement in CV systems in 
terms of the residual distributed entanglement. In this respect, 
the intriguing question arises whether a stronger monogamy 
constraint exists on the distribution of entanglement in many- 
body systems, which imposes a physical trade-off on the shar- 
ing of both bipartite and genuine multipartite quantum cor- 
relations. It would be important to understand whether the 
inequality [Eq. ©] holds as well for discrete-variable quc/its 
(2 < d < 00), interpolating between qubits and CV systems. 
If this were the case, the (convex-roof extended) squared neg- 
ativity, which coincides with the tangle for arbitrary states of 
qubits and with the Gaussian tangle for Gaussian states of CV 
systems, would qualify as a universal bona fide, dimension- 
independent quantifier of entanglement sharing in all multi- 
partite quantum systems. In such context, a deeper investi- 
gation into the analogy between Gaussian states with finite 
squeezing and effective finite-dimensional systems, focused 
on the point of view of entanglement sharing, may be worthy. 
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Appendix. — In this Appendix we prove the inequality in the 
right-hand side of Eq. d20l . From Eqs. ( fT6b and ( TTTb . we have 
f ± -j] = -(ab-c 2 _)[(ab-c 2 + )(ab-c 2 _)+l-a 2 -b 2 +2c + c-]. Using 
Eqs. (TBUIBl we find £ 2 > 77 so that ft = ht(n) = - ^/j? > 
(£_ > 0) and £ 2 < 77 so that f + < -^rj. Here, Eq. (fT3l implies 
£ + > c + (ab-c 2 _-\) > 0. However, £ + = implies det y Al A, = 
(ab - c 2 )(ab - c 2 ) = 1, which means that y Ai:Al is the CM of a 



4 



pure Gaussian state. This can be ruled out from the beginnign 
and therefore we have £ + > 0. Substituting y[rj > £ + into 
= £_ — -\Jrj, we obtain 

0<fi = -2c_ < 2 V 3 ^ - ^. (22) 

Next, from Eqs. ([Oi l and ( fT~8T >, we observe f + (a 2 + fr 2 )(c + + 
c-)(ab - c 2 - 1) = -(a - fV) 2 [c + - c-(ab - c 2 )] < to obtain 
£ < -(a 2 + b 2 ){c+ + c-)(ab - c 2 - 1) < 0. The last inequality 
is again due to Eq. dl3) . Hence, we obtain 

£, = h 2 (n) = 2(ab - c 2 )(a 2 + b 2 + 2c + c_) + f/ V?7 
> 2(a 2 + Z? 2 + 2c + c_) + f = £ . (23) 

Here, we have used Eq. ( fT3l , the inequality > > 0, and 
a 2 + b 2 + 2c+c- > 2(ab - c 2 ) > 0. Now, we observe that: 
- 4) = 2(a 2 + b 2 + 2c+c_)£ + + £ - 4£+ = -4c_ + 3(a 2 + 
b 2 )c--2a 2 b 2 c-+2abcl-2abc + + (a 2 +b 2 )c + (ab-c 2 ) + 8c 2 c + + 
4abc-c\-4c i _c\ > -4c -[a 2 + b 2 +2c+c--(ab-c 2 + )(ab-c 2 )] + 
3(a 2 + b 2 )c- - 2a 2 b 2 c- + 2abct - 2abc+ + (a 2 + b 2 )c+{ab - c 2 ) + 
8c 2 c+ + 4abc-c\ - 4c\c\ = -(a 2 + b 2 )c- + 2a 2 b 2 c- - 2abcl - 
2abc+ + {a 2 + b 2 )c+{ab - c 2 ) > -2abc- + 2a 2 b 2 c- - 2abcl - 
2abc+ + 2abc+{ab - c 2 ) = 2ab(c + + cJ)(ab - c 2 - 1) > 0, 
where we have used Eq. ( TT~4-b . Noting that > 0, we obtain 
C 2 > 4. Finally, Eqs. <E2) and (|23), with ^ > 4, yield: ^ 2 /^ 2 < 
C 2 /C 2 = -4c + c-/f 2 < -c + c_ . " " ■ 
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